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ABSTRACT 


Using a new color-color diagnostic diagram in the mid infrared built from 
WISE data, the MIRDD, we compare narrow emission-line galaxies (NELGs) 
that exhibit different activity types (star-forming galaxies, SFGs, and AGNs, 
i.e.,LINERs, Sy2s and TOs), as determined using one standard diagnostic di¬ 
agram in the optical (BPT-VO), with broad-line AGNs (QSOs and Syls) and 
BL Lac objects at low redshift (z < 0.25). We show that the BL Lac objects 
occupy in the MIRDD the same region as the LINERs, whereas the QSOs and 
Syls occupy an intermediate region, between the LINERs and the Sy2s.In the 
MIRDD these galaxies trace a sequence that can be reproduced by a power law, 
Fiy = z/“, where the spectral index, a, varies from 0 to —2, which is similar 
to what is observed in the optical-ultraviolet part of the spectra of AGNs with 
different luminosities. 

For the NELGs with different activity types, we perform a stellar population 
synthesis analysis, conhrming that their specihc positions in the MIRD depend 
on their star formation histories (SFH), and demonstrating that the W2 — W3 
color is tightly correlated with the level of star formation in their host galaxies. 
In good agreement with the SFH analysis, a comparison of their MIR colors 
with the colors yielded by energy distributions (SEDs) of galaxies with different 
activity types, shows that the SED of the LINERs is similar to the SEDs of the 
QSOs and Syls, consistent with AGN galaxies with mild star formation, whereas 
the SEDs of the Sy2s and TOs are consistent with AGN galaxies with strong 
star formation components. For the BL Lac objects, we show that their blue 
MIR colors can only be htted with a SED that has no star formation component, 
consistent with AGNs in elliptical-type galaxies. 

From their similarities in MIR colors and SEDs, we infer that, in the nearby 
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universe, the level of star formation activity most probably increases in the host 
galaxies of emission-line galaxies with different activity types along the sequence 
BL Lac^LINER^QSO/Syl^Sy2^TO^SFG. 


Subject headings: galaxies: active — (galaxies:) quasars: general — galaxies: Seyfert 
— galaxies: star formation 
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Introduction 


A large frac 


;ion of galaxies in t 


le nearby universe s 


low emission lines in their nucleus 


(Heel 

cman 

1980; 

Kennicutt 

1992; 

Ho 

1996; 

Ho et ah 

1997; 

York et ah 

2000; 

Stouehton et ah 

2002; 

Torres-Panaoui et ah 

2011 

). However, only a few percent of these galaxies exhibit 


broad lines similar to quasars, where we are sure that the nuclear activity is due to the 


accretion of gas on t o a supermassive black hole (SMB 


Lvnden-Bell 

1969; 

Soltan 

1982; 


Rawlings fc Saunders 


) at the center of their host galaxies 


199l|). 


In the 1980s, different studies demonstrated that by comparing the ratios of different 
spectral emission lines in the optical it is possible to devise standard diagnostic diagrams 
that allow to identify the dominant source of ionization of the gas in nar r ow em ission-line 


galaxies (NELGs). In their seminal paper. 


Baldwin. Phillips fc Terlevichl ( 1981 


BPT) only 


recognised three possible sources: massive young stars in Star Forming Galaxies (SFGs), 
the accretion of gas onto a SMBH in Active Galactic Nuclei fAGNs), and shock-wave 


heating due to supernovae winds or their remnants. Later on. 


Veilleux fc OsterbrockI (119871. 


VO) argued based on empirical evidence that only the two hrst mechanisms proposed by 
BPT are most likely to explain the line ratios in NELGs, and proposed simple criteria 
to distinguish between SFGs and AGNs in standard diagnostic diagrams. In particular, 
one diagnostic diagram that compares the line ratio [NII]A6583 /IIq; with the line ratio 
[OIII]A5007/II/5—hereafter, the B PT-VO diagram—pr oved to be one of the most useful 


tool to make this distinction (e.g. 


Stasihska et ah 


209611 . 


Thanks to l arge spectroscopic sur veys like the Sloan Digital Sky Survey (SDSS; 


York et af 


2000 : 


Abazaiian et ah 


20091) . we have now the possibility to apply the BPT-VO 


diagram to very large samples of NELGs at low redshift. What i s observed, then, is a 


continuous sequence t 


Kauffmann et ah 


2003 


r at takes the shape o: 


) the gr eek symbol v (IKewlev et ah 


Torres-Papaqui et al 


2001 


2012al ). the SFGs forming the left branch and 
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the AGNs the right one (cf., Fig. [T]). 


The nature of t he sequence produced by t 


well understood ("see 


Coziol 


1996 


Coziol et ah 


le SF Gs in the BPT-VO diagram is now 
, and referenc es therein). T 


1999 


due to a variation in spiral galaxies of the meta' 


Dopita & Eyans 

1986; 

Torres-Papaaui et al. 

2011 


l icity o f the gas flMcCall et al 


ms IS 


1985 


2012b ): the ratio [NII]A6583/Ha rises 


with the metallicity, while the ratio [OIII]A5007/H/1 decreases as the level of excitation of 


the gas diminishes with the temperat ure. However, the explanation 


is more ambiguous and controversial f|Viegas fc de Gonveia dal Pino 


2000 


Nagao. Maiolino fc Marconi 


or the AGN branch 


1992 : 


Gooke et al. 


20061) . The most natural assumption is that the high 


[NII]A6583/Hq; ratio seen in AGNs—their distinctive trait in the BPT-VO diagram—reflects 
a high abundance of chemical elements. But this would imply that the gas metallicity in 
AGNs is much higher than in SFGs, ignoring the cooling effect of metals in interstellar gas: 
as the metallicity grows, the gas temperature decreases due to cooling by emission, reducing 


the ratio [OIII]A5007/H/1, and not increasing it as it is observed o n the right 


of the i/-shape distribution formed by AGNs (see the discussion in 


Goziol et al 


branch 


20111 ). 


their metallicity ( 

Osterbrock 

1970; 

Storchi-^ 

lergmann & Pastoriza 

1989; 

Store 

li-Bergmann 

I99I; 

Hamann & Ferland 

199f 

1: Goziol et al. 

1999 , 

2011; 

Torres-Papaqni et al. 

2011, 

2012a) 


But this would imply that the formation of SMBHs at the center of galaxies and the 


at the same time 

Doknehaev 

1991 

Hamann V Ferland 

1992; 

Bovle & Terlevich 

1998; 

Franceschini et al. 

1999; 

Haiman et 

al. 

2004; 

Menanteau et al. 

2005; 

Lon & Jiang 

2008 


Haiman et al. 2007; Li et a 


2007 


Bennert et al, 


20081: Fanidakis et al. 2011 


201 3 I: 


: iTreister et al 


2011 


Torres-Papaqni et al 


2012a 


Floyd et al 


Torres-Papaqni et al 


20131 ) 


Another difficulty in explaining the AGN branch in the BPT-VO diagram is related 

























































































































































6 


to the Low ionization Nnclear Emission-line Regions (LINERs) and the Transition-type 


(Heel 

man 

( 

CO 

00 

2 _ ^ 

Ho et al. 

199c 

1; Ho 

1996; 

Kewlev et al. 

2001; 

Ho. FilioDenko & Sargent 

2003; 

Kauffmann et ah 

2003 

; ] 

kauffmann 

2009 

). That LINERs are genuine AGNs is 


snggested by the extreme position of these galaxies in the BPT-VO diagram: together 
with Seyfert 2 galaxies (Sy2s), LINERs form a continnous seqnence at the right end of 
the z/-shape branch (c.L, Fig. [1]), implying, according to photoionization models, that the 
ionizing sonrce in a L 


flOsterbrock fc Dahari 


INER 

is the same as in a Sv2, but intrinsical 

V weaker in luminosity 

1983: 

Wil 

Iner et al. 

1985; 

Goodric 

1 & 

keel 

1986; 

Ho et a . 

1993 



Schmitt et ah 


Starling et ah 


1997; 


2005 


Barth et ah 1998: Larkin et ah 


Kewlev et ah 


2006 


Maoz 


2007 


Masegosa et ah 


1998; Ho 1999; Sarzi et ah 2005 


201lh . However, and 


despite the BPT-VO diagram, many researchers still donbt the AGN natn r e of LINERs 

3 


snggesting, instead, that they are ei t her ionized by shoc ks flHeckma 


Veilleux et a 

1995 

Kim et al. 

1998; 

Veilleux et al. 

1999 

), 

1992 

Dopita 

1995; 

Maoz et al. 

1998 

). Some authors have 


1986 


Forbes 


1993 


19991) . or by very massive stars fIShields 


LINERs is ionized by stars mnch less massive than O and B stars, bnt mnch hotter: the 
so-called post-asymptotic giant branch stars (PAGBs). This appellation, nnfortnnately, 
is being used vaguely in the literature, as it could refer to pre-planet ary nebulae (the 


1971; 

Taniguchi et a 

20G 

0 

Terashima et al. 

2000; 

Flohic et al. 

2006; 

Stasihska e 

al. 

2008 

Gid Fernandes et al. 

2011 


Giraud et al. 

2011 

Singh et al. 

2013 

Panaderos et al. 

201c 

i). 


Determining the nature of TOs based on the BPT-VO diagram is also problematic. 


because s 


ar formation is very inten se in these gal axies, while it i s mu ch lower in the 


LINERs flTorres-Papaqni et ah 


2013 ). This has led 


Kewlev et ah 


(1200 ll ) to propose new 


boundaries in the BPT-VO diagram between AGNs and SFGs, suggesting that TOs form a 
buffer zone (cf.. Fig. [1]), where star formation is still the predominant source of ionization 
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of the gas. Consequently, the most accepted interpretation for TOs seems to be that they 


are AG Ns hidden beh i nd ex t ended regions of intense star formation 


happening near 


heir 


center ( 

Ho et al. 

1993, 

1997; 

Veron. Goncalves & Veron-Gettv 

1997; 

Goncalves et al. 

1999 

Ho. FilioDenko & Sargent 

2993; 

Torres-Panaoui et al. 

2912a 

^ _ 2 

913 

)• 


On the other hand, explaining why the TOs differ from the LINERs and the Sy2s 
is not obvious, because all these galaxies trace an apparently continuous sequence in the 
BPT-VO diagram. One possible explanation is that TOs are AGNs at an earlier stage 
of evolution, whence the name “transition-type objects” , since they would be expected 


to eventually transform into eithe r Sy2s or LINERs when t 


le level o 


would have decreased sufficiently 


Vila-Vilaro 


2000 


Sarzi et ah 


2005 


Osterbroc 


c fc Shaw 


Lee et ah 


2007; 


1988; 


Wu et a. 


Haan et ah 


2008 


star 


1998 


brmation 


Maoz 


Wang k. Wei 


1999 


2008 


2999I: 

Ghen et al. 

2919; 

Garnineti et al. 

2912; 

Torres-Panaoui et al. 


Montero-Dorta et al 


2013|). 


The ambiguities as to the nature of LINERs and TOs, as inferred from the BPT-VO 
diagram, is a serious drawback in our understanding of the AGN phenomenon. It suggests 
that standard diagnostic diagrams based on emission-line ratios could be misleading, 
and that conclusions on the nature of AGNs that were classified using this method are 
possibly erroneous. Moreover, because we cannot apply the BPT-VO diagram to BLAGNs, 
establishing a connection between the different AGN types observed in the nearby universe 
is a complicated issue. What is needed, therefore, are alternative means by which we 
can verify independently the classification of the activity of NELGs as obtained using the 
BPT-VO diagram, and that could be used to compare NELGs with BLAGNs directly, i.e., 
by comparing the same character istics. This can no w be done thanks to the Wide-held 


Infrared Survey Explorer (WISE; 


Wright et al. 


291911 












































































































WISE is an all-sky photometric survey in four wave bands in the mid infrared (MIR): 
3.4 fim (Wl), 4.6 fim (W2), 12 fim (W3) and 22 /im (W4). The WISE All-Sky Release 


Source Catalog, which we use in this study, contains astro metry and photo metry for 


563,921,584 objects detected on the Atlas Intensity Images (jCutri et ah 


201311 . This survey 


has already produced important results, showing, in particular, that e mission-line galaxie s 


with different activity types have distinctive MIR characteristics (e.g. 


Stern et ah 


2013 


2012 


Yan et ah 


Ma 


eos et ah 


2012 


Donoso et ah 


2012 


Jarrett et al. 


Clemens et al 


201311 . Following the same line of research. 


20 


Coziol et al 


1 


2011 


Assef et al. 


(1201411 have 


recently devised, using WISE data, a new diagnostic diagram, the MIRDD, that allows to 
separate empirically NELGs with different activity types, in a way that is fully consistent 
with the classification obtained from the BPT-VO diagram. The MIRDD consists in 
comparing the W2 — W3 colors with the W3 — W4 colors, and although the separation 
criteria are empirical, it is suspected, as it is also the case for the BPT-BVO diagram, that 
the reason why such segregation is possible in the MIR is because these two colors vary 
systematically as the contributions of the two main components forming the spectral energy 
distribution (SED) in the MIR, an AGN and star formation in the host galaxies, change 
with the activity type. 

For example, in the MIRDD the AGN nature of LINERs seems now much more 
compelling. This is because their MIR colors are signihcantly different from those of SFGs 
and PAGBs, implying that although massive stars or white dwarfs might still be assumed 
to ionize the gas in LINERs, these stars obviously cannot explain their MIR emission, 
whereas an AGN can explain both. Moreover, the LINERs and Sy2s trace in the MIRDD 
a continuous sequence that is consistent with the color sequence produced by power laws 


with different spectra 


Osterbrock fc Ferland 


indexes, in 


20061; 


Maoz 


good agreement with photoionization models (e.g.. 


20081 ). whereas SN shock-wave is not a mechanism that is 


capable of heating dust. Also, the W2 — W3 colors of the NELGs in the MIRD were found 
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to become gradually redder along the sequence LINER—)-Sy2—)-TO—)-SFG, which suggests 


that this color is mostly sensitive to the presen' 


galaxies of emission-line galaxies flDonoso et ah 


leve 


2012 


of star formation activity in the h ost 


Torres-Papaqui et al. 


2012 a, 


201811 ■ 


Consequently, we here propose a new study that should allow us to better understand 
how the MIRDD works, and gain new information on how BLAGNs are connected with 
AGNs in NELGs. This article is organized in the following way. We introduce our sample 
of nearby {z < 0.25) NELGs and BLAGNs in Section [2], and use the BPT-VO diagram 
to separate the NELGs into six different activity types. In Section 13. 1[ we discuss how 
these galaxies are differentiated in the MIRDD, and support this discussion by doing a 
statistical analysis of the MIR colors. In Section 13.1.11 we examine what are the effects on 
our statistical analysis of having different distributions in redshift and different number 


of resolved galaxies in our sample. In Section 13.21 we perform a new stellar populatio n 


synthesis analysis on the I 

'lELGs (e.g 

Goziol et ah 

2011 ; 

Torres-Panaani et al. 

2012 a 

Plauchu-Fravn et al. 

2012 ; 

Torres-Panaani et al. 

2013 

), deducing their star formation 


histories (SFHs), to check how well this characteristic can explain the positions of galaxies 
with different activity types in the MIRDD. We then proceed by studying the relation 
in NELGs between the W2 — W3 color and the intensity of star formation in their host 
galaxies. In Section 13.31 we complement the SFH analysis by studying how the changes of 
colors in the MIRDD are related to the different components that are forming the SEDs of 
emission-line galaxies with different activity types. Finally, in Section l3^ we compare the 
MIRDD of BLAGNs and BL Lac objects with the MIRDD of NELGs, and in Section 13751 
we compare their SEDs. In our discussion. Section IH we consider the consequences of 
our comparison of BLAGNs with NELGs in the MIR, and summarize our results and 
conclusions in Section |5l 
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2. Samples and MIR data 


The NELGs used to compare the BPT-VO diagram with the MIRDD were selected 
from the main spectroscopic catalog of the Sloan Digital Sky Survey Data Release 7 (SDSS 


DR7; 


Abazaiian et al. 


20091) . which represents the completion of the SDSS-II project. The 


SDSS DR7 catalog includes the spectra of 9.3 x 10^ galaxies, 1.2 x 10^ quasars and 4.6 x 10® 


rjil. 


starqj. Only the objects identihed as galaxies were considered for our study, keeping all 


those that have a redshift ^ < 0.25. After correcting for t 


he red 


templates produced by STARLIGHT fiGid Fernandes et al. 


shifts, stellar population 


20051) were subtracted from 


the data, leaving pure emission spectra. For our study, we selected only the NELGs that 


have optical emission line ratios, [NII]A6583/Ha and [OIII 


ratio S/N > 3, and S/N > 1 0 in the adjacent continuum flRrinchmann et aLll2004l l. which 


A5007/H/3, with a 


signa l-to-noise 


amounts to 229,618 NELGs flTorres-Papaqui et al. 


2012 


To retrieve the MIR data, we cross-correlated the positions of the NELGs, as found in 
SDSS DR7, with the positions of the objects in the WI SE All-Sky Release Source Gata log. 


This was done using the X-Match pipeline in VizieR flOchsenbein. Bauer &: Marcout 


2000 lj. applying a se arch radius of hve arcseconds around the position of each galaxy (e.g.. 


Glemens et al. 


20131) . For each match, the algorithm determines the angular distance from 


the candidate to the target. By plotting this angular distance, we found distributions that 
peak around 0.2 arcsecond, with relatively very few galaxies (~ 11%) with a distance larger 
than 1 arsecond. By reviewing the WISE identihcation, we eliminated all the duplicate 
matches (which represent much less that 1% of the candidates), keeping only those that 
have the smallest angular differences. Then, selecting only the NELGs that have WISE 
fluxes with S/N > 3 (corresponding to quality flags, ph_qual^ equal to A or B) in all of 
the four WISE MIR bands, and that have a contamination and confusion flag that is clear 


^http: / / WWW. sdss.org/dr7 / start / aboutdr7.html 
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[cc-flags = 0), we obtained 133,095 NELGs (i.e., 58% of our original spectroscopic sample). 


Using the BPT-VO diagram, we classified our sample of MIR NELGs according to 
their activ ity type. The results a re sho wn in Fig. [1] The different activity types are those 


dehned in 


Torres-Papaqui et ah 


(120131) . The NELGs are separated in three main types. 


the SFGs (66.7% of the MIR NELGs), the TOs (22.7%), the Sy2s (7.0%) and the LINERs 


(3.6% ). Our MIR sample of nearby NELGs cle arly favors ga^ 


is relatively high (similar to what was found by 


Donoso et ah 


axies where star formation 


2 OI 2 I) . We also separated 


the SFGs and TOs in two subtypes. Gompared to the SFGls (56% of the SFGs), the 


SFG2s are star 
explanations in 


brming galaxies that show a n excess of nitrogen due to starburst winds (see 


Torres-Papaqui et ah 


2012bf) . whereas the T02s (74% of the TOs) differ 


from the TO Is because of their lower leve 


LINERs than Sy2s flTorres-Papagni et ah 


of ex citation, which makes them more similar to 


2013 ). Our hnal WISE selected, NELGs sample 


is thus composed of 49,724 SFGls, 39,026 SFG2s, 7,752 TOls, 22,671 T02s, 9,361Sy2s and 
4,741 LINERs. 

The BLAGNs in our sample come from the 1 3th e dition of the catalog of quasars and 


active nuclei published by 


Veron-Gettv fc Veronl fl2010l) . which is also accessible through 


VizieR. From this catalog, we selected all the entries identihed as QSOs or Syls (more 
specihcally for the Syls, only those identihed as SI in the catalog). Although it has become 
customary to separate quasars into type 1 and type 2, we have made no such distinction 


in ou r study, because this would have implied doing a followup in X-rays (e.g. 


Derry et ah 


20031 ). which is beyond the scope of the present study. Also, although BL Lac objects 


are not BLAGNs (since their spectra show no emission lines), we have included them in 
our study because, due to their non-thermal continua in the optical, they are presumed 
to represent the “purest” AGN cases p ossible, and, consequently, they can be i mportant 


points of reference in the MIRDD (e.g.. 


Plotkin et al. 


2012 


Massaro et al 


20131) 
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Like for the NELGs, we restricted the redshift of our candidates to ^ < 0.25, obtaining 
three initial spectroscopic samples of 7,038 Syls, 3,510 QSOs, and 169 BL Lac objects. The 
MIR data from the WISE All-Sky Release Source Catalog were retrieved in the same way 
as for the NELGs. Comparing the angular distances we found wider distributions than for 
the NELGs, with a maximum near 1 arcsecond (compared to 0.2 arcsecond for the NELGs), 
but with only a few percent of candidates with angular distances larger than 2 arcseconds. 
After eliminating all the duplicated matches, and applying the same selection criteria as 
the NELGs for the MIR fluxes, we ended up with 5,077 Syls, 2,323 QSOs, and 55 BL Lac, 
which represent, respectively, 72%, 66% and 33% of the initial spectroscopic samples. 

In Fig.|2]we show the normalized distributions in redshift of all the selected emission-line 
galaxies, as separated by activity types. As it should have been expected, the NELGs are 
better represented at low redshifts than the BLAGNs. The distributions of the SFGs and 
TOs are similar, showing a maximum at z ^ 0.07 and a long exponentially decreasing tail, 
that stays well populated up to 2 ; ~ 0.14 for the two SFGs and T02s, but decreasing less 
rapidly with the redshift in the case of the TOls. The SFG2s also differ from the SFGls, 
presenting a higher fraction of galaxies at low redshifts. The distribution for the Sy2s is 
similar to the distributions of the TOls. The distribution for the LINERs is similar to 
the distribution of the SFG2s, showing an almost constant fraction of galaxies from z = 0 
to z = 0.07, the fraction of galaxies decreasing rapidly at higher redshift. For the BQSOs 
we see the inverse trend, their fraction increasing with the redshift. The QSO distribution 
reaches its maximum only at z = 0.25, while the distribution of the Syls reaches its 
maximum at much lower redshift, their fractions remaining constant over the range ^ = 0.08 
to z = 0.25. The difference between the Syls and QSOs is as expected from the classihcation 
criterion, which is based on a difference in absolute magnitude: consequently, the fraction 
of BLAGNs classihed as QSOs instead of Syls is naturally expected to increase with the 
redshift. The BL Lac objects show no obvious trend with the redshift. But, this is most 
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probably due to small statistics (there is only 55 objects in this sample). 

From the different redshift distributions it is expected that a greater number of NELGs 
will be spatially resolved in WISE. Indeed, revising the values of the flag exLflags, we find 
that 66% of the LINERs are resolved, compared to 28%, 35% and 23% of the Sy2s, TOs 
and SFGs, respectively. Similarly, we find that a high fraction of BL Lac objects are also 
resolved, ~ 47%, which is relatively high considering their small number and homogeneous 
distribution in redshift. In comparison, only 17% of the Syls and 9% of the QSOs are 
resolved. In the next section, we will examine in more detail the possible consequences of 
these differences on our analysis. 


3. Results 


3.1. Tracing the MIRDD for nearby NELGs with different activity types 


We show the MIRDD for the NELGs in Fig. [3] for the Sy2s and the LINERs, in 
Fig m for the TOS and in Fig. [5|for th e SFG s. The different regions of the MIRDD were 


determined empirically in 


Goziol et ahl (120141) . The analytical relations identifying these 


different regions are the following: ii x = W2 — W3, and y = W3 — W4, then the separation 
between AGNs and SFGs is given by the expression y = 1.6a; — 3.2, and the separation 
between high and low star formation activity is given by the expression y = —2.0a; + 8.0. 
One can see that on the AGN side we find mostly Sy2s and LINERs, whereas on the side of 
low star formation activity the majority of the NELGs are LINERs. 


The MIRDD and the BPT-VO diagrams yield highly consistent results for the LINERs 
and Sy2s. In both diagrams these two AGNs trace a continuous sequence, the Sy2s differing 
from the LINERs in the optical by their higher level of ionization, and, consequently in 
the MIR, by their slightly higher dust temperatures. Moreover, we find that the colors 
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of the LINERs and Sy2s are consistent with the colors produced by different power laws, 
Fy = I/", the spectral index, a, varying between 0 and —1.5 in the LINERs and between 
— 1 and —2 in the Sy2s. As they are, these different values are consistent with what is 


obse rved in the optical-ultraviolet region of 


le.g. 


Osterbrock fc Eerland 


2006 


Maoz 


20081 ) 


he spectra of AGNs with different luminosity 


Note that the position of the LINE Rs in the MIRDD is different from the position 


of the LINERs as found in Figure 12 in 


Wright et al 


120101 ) . which has been taken as 


reference by m any authors in the li terature. The reason for such a difference is because the 


LINERs in the 


Wright et al. 


1120101) sample are ultra-luminous infrared galaxies ( ULIRGs), 


which are quite different from the “standard” 


1119801) . This difference is well documented in 


INER s as de hned in the optical by 


Heckman 


Gozioll ljl996[l . where it is shown that what 


distinguishes these two types of LINERs is their differ ent levels of s 


as one can verify, the ULIRGs classihed as LINERs in 


Wright et ah 


; ar form ation. Indeed, 
1120101 ) occupy in the 


MIRDD the same position as the TOs and SFGs in Fig. |5l It is important to stress this 
difference, because it demonstrates that although emission-line galaxies may have similar 
spectral characteristics in the optical, their MIR colors, and their W2 — W3 colors in 
particular, may differ signihcantly, depending on the different levels of star formation in 
their host galaxies. 

Similarly, the different positions occupied by the Sy2s and LINERs i n the MIRDD ar e 


consistent with a higher level of star formation in Sv2s t 

han in LINERs ( 

Lawrence et al. 

1985; 

Dultzin-Hacvan et ah 

1988; Torres-Panaani et al. 

2013: Rosario et al. 

2013 


Goziol et al. 

2014 

). This may also explain their different power laws, the steeper spectral 


i ndex for the Sy2s implying a n excess in IR emission, which could be due to star formation 


llTorres-Papaqui et al. 


2 m2bh . 


The MIRDDs for the TOs in Fig 0] is also in good agreement with the BPT-VO 
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diagram. In both 


SFGs and AGNs f Kewlev et al. 


diagrams, the TOs seem to occupy a sort of buffer zone between the 


200ll) . The fact that this is due to star formation is more 


obvious in the MIRDD, where the W2 — W3 colors for the TOs are significantly distinct 
from those of the LINERs (contrary to the continuous sequence shown in the BPT-VO 


diagram). The MI RDD of the TOs clearly emphasizes the high leve 


their host galaxies (IKewlev et al. 


2001 


Torres-Papaqui et al. 


2012 al. 


s of star formation in 


2013h . 


The differences between the TOs and SFGs are better appreciated by comparing Fig. 0] 
with Fig. O The SFGls have slightly redder W2 — W3 colors and bluer W3 — W4 colors 
than the SFG2s. The highest similarities are between the TO Is and SFG2s, which, again, 
is in good agreement with the relative positions occupied by these galaxies in the BPT-VO 
diagram. As expected, the colors of the SFGs are clearly distinct from those produced by a 
power law. On the other hand, none fall on the black body curve. 


In general, we distinguish in the MIRDD a dehnite trend for the W2 — W3 color to 
become redder along the sequence LINER—)-Sy2—)-T02—)-T01—)-SFG2—>-SFGl. This trend 
is obvious in Fig. where we present the box-whisker plots for the MIR colors of the 
different NELGs in our sample. These plots include notches around the medians that are 
not overlapping, which suggests that the medians are significantly different. This is not the 
case for the W3 — W4 colors shown in Fig. [6b. The SFGls and SFG2s have similar medians, 
and so have the TOls and T02s. There is, consequently, a weaker trend in W3 — W4 
colors, the NELGs becoming bluer, by only 0.5 mag, along the sequence Sy2—)-TO—)-SFG. 
The LINERs do not follow this sequence, being slightly bluer than the SFGs. 


These differences in colors are confirmed statistically in Fig. [TJ where we present 


(Hothorn. Bretz & Westfall 

20 Q§; 

Herberich. Sikorski & Hothorn 

201 Q) 


20101 ). A more explicit 


explanation of this parametric statistical test, and how to interpret its results, can be found 
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in 


Torres-Papaqui et al.l (120131) . In summary, the test compares the mean MIR colors, 


family-wise, i.e., comparing simultaneously all pairs of subsamples possible. For example, 
in Fig. [7^ the test for the W2 — W3 colors conhrms that the LINERs are bluer by ~ 1.2 
mag, than the SFGls, while the smallest differences are found between the T02s and TOls 
(smaller than 0.1). In Fig. [7b, the test also conhrms that the differences between the mean 
colors are much smaller in W3 — W4 (huctuating around zero by ~ 0.3 mag) than in 
W2 - W3. 


3.1.1. Effect of the resolved galaxies on our statistical analysis 


Our statistical analysis clearly shows that the W2 — W3 color is a very sensitive 
parameter in the MIRDD, varying monotically by ~ 1.2 mag along the sequence 
LINER—)-Sy2—)-T02—)-T01—)-SFG2—)-SFGl. This suggests that the W2 — W3 color must 
be related to some intrinsic physical characteristic that is systematically varying with the 
activity type of the galaxies. On the other hand, we have also found that our samples 
contain different fractions of resolved galaxies at different redshifts. In particular, our 
sample of LINERs turned out to include a signihcantly larger number of resolved galaxies 
at lower redshifts than the other samples. Thus, it seems important to examine what is the 
importance of these differences on the results of our statistical analysis. 


According to the WISE explanatory supplemento, the photometry for resolved sources, 
as identihed by the flag ext-flags, underestimates their fluxes. Physically, this means that 
for the resolved galaxies only their inner regions are included in the photometric aperture. 
And even more so, if the galaxies are at low redshift. On the other hand, because the 
contamination and confusion flags are clear {cc_flags= 0), and the quality flags {ph_qual) are 


^http: / / wise2. ipac .caltech.edu/docs/release / allsky / expsup / 
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equal to A or B for all the resolved galaxies, we are sure that the data in these regions are 
of good quality. What is left to establish, therefore, is the effect on our statistical analysis 
of having only partial fluxes for the resolved galaxies. 


To determine this effect, we have compared the colors of the unresolved galaxies 
{ext_flags = 0) with those of the resolved galaxies [ext^flags ^ 0). What we found, is 
that, in general, i.e., independently of the activity type, underestimating the fluxes makes 
the MIR colors bluer. There are two possible causes for this “blueishness effect” (BEff). 
One is an observational bias: in the MIR, older stellar p opulations haye typ i cal blue colors 


fe.g., the colors of a K-ty pe star in the MIRDD; see also 


2012 


Jarrett et ah 


Wright et ah 


2010 


Plotkin et ah 


20131 1. and since stellar populations in the center of galaxies tend to 


be older, on average, than in its periphery, then, the resolved galaxies are expected to be 
bluer than the unresolved galaxies. The other cause is physical: bec ause galaxies with 


different morphologies do n ot form their stars in the same way (e.g.. 


Torres-Papaqui et al. 


Coziol et ah 


2011 


20131) . then galaxies with early-type morphologies will be expected to 


have bluer colors, on average, than galaxies with later-type morphologies. 


To distinguish between the two causes described above, we have first quantified the 
BEff comparing in each activity-type samples the mean colors of the unresolved galaxies 
with the means for the total samples. We found that the BEff for the W2 — W3 and 
W3 — W4 colors are of the order of 0.35 and 0.32 mag for the LINERs, compared to 0.04 and 
0.06 mag for the Sy2s, 0.05 and 0.07 mag for the TOls, 0.09 and 0.10 mag for the T02s, 
0.02 and 0.00 mag for the SFGls, and 0.04 and 0.03 mag for the SFG2s. For comparison 
sake, the BEff for the BL Lac objects are 0.16 and 0.01, and practically zero for the Syls 
and QSOs. So, only the LINERs (and BL Lac objects) are significantly affected. This 
implies that the BEff cannot be explained by an observational bias, but that LINERs must 
be physically different from the other galaxies. 
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To examine this point in more detail, we compare in Fig. |8]the distributions of the 
W2 — W3 colors, which is the parameter that shows the largest variations, as a function of 
the redshift. In each activity-type samples the galaxies were separated into four groups (0, 
1, [2,3,4] and 5), as given by the WISE photometry flag {ext_flags): where 0 indicates a 
point source, and the degree of resolution increases from 1 to 5. As expected, we see that 
the fraction of resolved galaxies increases at low redshift. However, the BEff becomes more 
important at low redshift only in the LINERs, which confirms that this effect is caused by 
a physical difference. 


To support this interpretation, we show in Fig. IHlfor the NELGs the distributions of 
their W2 — W3 colors as a function of their Petrosian radii in kpc, as measured in the SPSS 


z band (assuming a pure Hubble flow, with Hq = 75 km s ^ Mpc 


Torres-Papaqui et ah 


2012 aJ). It can be seen that the level of resolution increases with the physical size of the 


NELGs. In general, the LINERs contain more extended galaxies than the other NELGs, 
and the BEff systematically increases with their size. Moreover, for comparable ranges in 
redshift and comparable ranges in galaxy size the LINERs always show bluer colors. This 


can only be explainec 


in the other NELGs (ITorres-Papaqui et al. 


assuming, on average, older stef 


2012 a, 


ar populations in the LINERs than 


20131) . Gonsequently, this result supports 


our statistical analysis, suggesting that the differences in W2 — W3 colors between NELGs 
with different activity-types are related to different stellar populations in their host galaxies. 


Also according to the WISE explanatory supplement, it is suggested that in the case 
of resolved galaxies one should use the magnitudes obtained by integrating the fluxes using 
the photometry apertures for 2MASS (i.e., using wlGmag instead of wlmag). To verify 
how this would affect our statistical analysis, we have consequently retrieved the available 
aperture photometry measurements in the WISE All-Sky survey using the GATOR catalog 
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query service at IPAC/IRSAlf|. We have found data for 83% of the resolved LINERs, 67% of 
the Sy2s, 57% of the TOls, 66% of the T02s, 47% of the SFGls, sand 55% of the SFG2s. 
For comparison sake, we also retrieved the data for the BL Lac objects and BLAGNs, and 
found magnitudes for 67% of the resolved Syls, 59% of the QSOs and 31% of the BL Lac. 
Gomparing the colors, we have found that the trend for the resolved galaxies to become 
bluer signihcantly increases. This is shown in Fig. [101 where we compare the colors obtained 
using wlGmag (W2G — W3G) with the colors using wlmag. One can see also that this 
effect does not depend on the different photometry flag (we used wAgflag). 

Using the wlGmag^ the BEff for the W2G — W3G and W3G — W4G colors would be 
0.52 and 0.41 mag for the LINERs, 0.57 and 0.24 mag for the Sy2s, 0.62 and 0.24 mag for 
the TOls, 0.60 and 0.25 mag for the T02s, 0.58 and 0.19 mag for the SFGls, and 0.60 and 
0.25 mag for the SFG2s. Similarly, we found a BEff of 1.1 and 0.24 mag for the BL Lac, 
0.43 and 0.25 for the Syls, and 0.62 and 0.29 for the QSOs. It is interesting to note that 
the resolved BLAGNs show a similar effect as the NELGs, suggesting that even in luminous 
AGNs, the MIR colors may be sensitive to a variation in stellar populations. In particular, 
using the wlGmag we have found a few BL Lac objects to have negative colors, which is 


typical of inact ive, elliptical ga 


these galaxies (iKotilainen et al 


axies, which is as expected from what is kno wn a priori of 


1998 


Scarpa et al. 


2000 


Falomo et al. 


20i4 . 


In general, taking into account that the fraction of resolved galaxies varies signihcantly 
with the sample, we have calculated (using a weighted mean) that only the LINERs would 
be signihcantly ahected, becoming bluer, in both colors, by ~0.2 mag compared to the 
other galaxies. Gonsequently, using the wlGmag instead of wlmag would change nothing 
to our statistical analysis. 


^http://irsa.ipac. caltech.edu/about.html 
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3.2. Variation of star formation in NELGs with different activity types 


In the previous section, we have shown that the systematic variation of W2 — W3 
color of NELGs w ith different activity types must be related to a variation in their physical 
characteristics. In ICoziol et al.l fj2014i) it was suggested that this characteristic is the present 


level of star f ormation in the host galaxies of t 


le NELGs. This was also s uggested by other 


authors fe.g.. 

Mateos et al. 

2012; 

Jarrett et al. 

2013; 

Rosario et al. 

2013 

). In particular. 

Donoso et al. 

(2012 

) have shown that the W2 — W3 is a valuable indicator of star formation 


activity, when no other data are available. This suggests that the W2 — W3 color sequence, 
LINER—)-Sy2—)-T02—)-T01—)-SFG2—)-SFGl, traced by the NELGs with different activity 


the MIRDD, may be c 

ue to an increase in the level of star formation in their host 

(Torres-Papaaui et al. 

2012a. 

2013 

)• 


It is possible to verify directly this interpretation for the W2 — W3 color sequence 
by applying a stellar popul ation synthesis analysis usi ng STARLIGHT to the SDSS 


spectra of the NELGs (e.g.. 


Torres-Papaqui et al. 


2 OI 3 I) . The stellar population synthesis 


method consists in quantifying in the optical spectra of galaxies the contributions of 
stellar populations with different luminosities or masses. One can then transform these 
contributions in terms of star formation rates (SFRs), and obtain in this way a snapshot 
of how a galaxy formed its stars over its lifetime, which constitutes its star formation 


hi story ( SFH). Furthe r details on how the SF L 


m 


Goziol et al. 


bm, 


Torres-Papaqui et al. 


is ob 


f 2012al ) and 


ained with STARLIGHT can be found 


Plauchu-Fravn et al. 


( 1201211 . 


Note that because the SDSS spectra cover mostly the inner parts of the galaxies, this 


l egitimates a comparison of 


(iTorres-Papaqui et al. 


2m2bl l. 


he SFHs with the colors obtained using only wlmag in WISE 


In Fig. [m we show the SFHs of the NELGs in our sample. We hnd that the SFH varies 
systematically with the activity type: the LINERs and Sy2s show their highest SFRs in the 
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past, star formation in the TOs is almost as high today as in the past, while it is highest 
at the present time in the SFGs. To make these differences ev en clearer, we have ex t racte d 


from the SFHs two characteristic SFRs as dehned in Fig. 6 of 


Torres-Papaqui et al. 


(120131): 


SFRyoung corresponds to the maximum SFR of the most recent star formation event, and 
SFRoid corresponds to the maximum SFR in the past. The means for the characteristic 
SFRs in NELGs with different activity types are compiled in Table [U In Fig. [T^] we 
compare their box-whisker plots. In Fig. \\2k we hnd that SFRyoung increases along the 
sequence LINER—)-Sy2—)-T02—)-T01—)-SFGl—)-SFG2. Note that SFRyoung is higher in the 
SFG2s than in the SFGls, w hich is explained by recent bursts of star formation in SFG2s 


(iTorres-Papaqui et al. 


2012 bl ). In Fig. [T2 ]d we find that SFRom decreases monotonically 


along the sequence LINER—)-Sy2—)-T02—)-T01^SFG2—)-SFGl. The differences between the 
SFRs as observed in the box-whisker plots are all confirmed statistically in Fig. [T31 showing 
the conhdence intervals for the max-t tests. 


In Table [T] we also give the age of the stellar populations producing the maximum SFR 
in the past, which is consistent with the time in the past when the star formation activity 
was at its maximum, tsFMAx- In Fig. [HI we find that tsFMAx increases in parallel with 
SFRoid- From this figure we also infer that, although SFRqw are comparable in the LINERs 
and the Sy2s, the former formed the bulk of their stars earlier than the latter. Similarly 
in the SFG2s, despite their recent burst of star formation, the star formation in the past 
reaches a maximum long before than in the SFGls, implying that the stellar populations in 
the SFG2s must be older, on average, than in the SFGls. 

In Fig. [island Fig. [16] we now compare the two characteristics SFRs with the two MIR 
colors used in the MIRDD. In Fig. [15^, except for the SFGls and SFG2s as noted before, 
we hnd that the W2 — W3 color is indeed well correlated with SFRyoung- The W2 — W3 
color of NELGs are getting redder by ~ 1.2 mag as the present level of star formation 
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increases in their host galaxies. In Fig. [16^, we see the inverse sequence, W2 — W3 getting 
redder as SFRom increases. On the other hand, in Fig. IT^ and in Fig. [TBb . we hnd no such 
correlations: except for the Sy2s, the W3 — W4 of NELGs are very similar. 


The above comparisons conhrm the high sensitivity of the W2 — W3 color to the SFH 
of the NELGs. The position of a NELG in the MIRDD depends not only on its nature as 
AGN or SFG, but also on its SFH, and, consequently, on the present level of star formation 
in its host galaxy. 


These results suggest that the main reason why the MIRDD (and BPT-VO diagram) 
allows to separate NELGs with different activity types is because they h ave different SFH s, 


i.e., galaxies with c 

ifferei 

it activitv tvpes formed 

dieir stars differentIv ( 

Dc 

kuchaev 

1991; 

Hamann & Ferland 

1992 

Bovle & Terlevich 

1998; 

Franceschini et al. 

1999; 

Haiman 

et al. 


2004: Menanteau et ah 2005: Lou fc Jiang 


2008: Haim an et a. 


2007: Li et al. 


20071: 


Bennert et al. 20081 


Torres- Papaqni et al 


anida 


2012 a 


ds et al. 2011: Treis 


Flovd et al. 


2013 


er et al. 2011: Donoso et al. 2012 


Torres-Papaqni et al. 


20131). 


3.3. SED study for the NELGs 


In the previous section, we have demonstrated that the W2 — W3 colors of the NELGs 
vary systematically with the SFHs of their host galaxies. However, the MIRDD also 
separ ates AGNs from S FGs, which suggests that the AGN component may also play some 


role (IGoziol et ah 


20141) . To determine what is the importance of the AGN component 


in the SEDs of emission-line galaxies with different activity types, we have compared the 
W2 — W3 colors of the NELGs with the colors produced by MIR SED templates as found 
in the literature. This study is important, not only because it should allow us to verify 
independently our interpretation based on the SFHs, but also because the same analysis 
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can be applied to BLAGNs, for which a stellar population synthesis analysis is obviously 
not possible. 

For our study we choose six represe ntative SEDs as show n in Fig. [T71 The hrst hve 


SEDs come from the library provided by 


Polletta et ah 


11200711 . The SEDs #1 and #2 were 


chosen to represent SFGs. These are two typical starburst galaxies, M82 and Arp 220, that 
due to their high level of star formation are also ultraluminous in infrared. The SEDs ^3 
and 7^4 represent, respectivley, a typic al Sy2 and Syl (Mrk 231). The SED #5 was chosen 


among the different SEDs proposed by 


Polletta et ah 


(120071 ) to represent a 


Our t emplates sample in Fig. [T7]is complemented with the SED ^6, from 


typical QSO. 


Leipski et ah 


(1201411 ■ representing a QSO with a strong emission in the FIR. 


For normal galaxies 


SED is straightfo r ward (iKennicutt et al 


(with no dominant AGN component) 


Sheth et al. 

2010; 

Jarrett et al. 

2 OI 3 I 


2003 


Murakami et al. 


;he in 

;erpretation 0 

■ the 

2007; 

Wrig-ht et al. 

2010 


201311 . Govering the wavelengths from ~ 1/rm to lOO/im, 


we can distinguish three different regions. Evolved stellar populations, dominated by low 
mass stars, produce the light at wavelengths lower than A ~ 3/im, forming a bump in 
emission that culminates at about l/rm. At longer wavelengths than A ~ lO/im, there is 
another emission bump that reaches a maximum at about lOO/im, and which is due to 
hot dust heated by young, massive stars in active regions of star formation. Finally, in 
between these two regions, 3 < A < 10/rm, there is a sort of valley, that is also very rich 
in spectral structures, showing either unidentihed IR bands (UIRBs) in normal galaxies, 


or PAH emission due to 

lot stars in starburst 

0 

0 

0 

CM 

Galzetti et al. 

2007; 

Hanami et ahl 

2012) 


;-like galaxies (jXu et al. 


1998; 


Laurent et al. 


201211 . Gonsequently, the infrared SED of normal 


galaxies is sensitive to stellar light from the evolved stellar populati ons, as well as t o 


low-temperature processes from t 

le ISM 

, and star formation activity ( 

Kennicntt et al.l 2003: 

Galzetti et al. 

2007; 

Wright et al. 

201 J; 

Sheth et ah 

2010 

Mateos et ah 

2012; 

Donoso et ah 
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2012 ; 

Rosario et al. 

2013; 

Jarrett et al. 

2013) 


These different components of the MIR SEDs of galaxies, as described above, can 
easily be distinguished in Fig. [T7]for emission-line galaxies with different activity types. For 
example, in SEDs ^ 1 and ^^2 we clearly see the two bumps at short and long wavelengths 
and the valley in between. Note that the depth of the valley is more or less pronounced, 
depending on the strength of the PAH emission (obvious in M82, but also in SED ^ 3 for 
the Sy2) and the presence of UIRBs, which are common features produced by different 
stellar populations in star forming galaxies. On the other hand, in the SED of AGNs, like 
SED #3 (Sy2) and SED #4 (Syl), the MIR valley more or less disappears. This is due to a 
high temperature dust component, which is growing in intensity between 1 /im and 10 /rm. 


and is interp reted to be re 


SMBH (e.g.. 


Laurent et al 


ated to dusty structures, near or at their center, 


2000 


Hanami et al. 


2012 


Roseboom et al 


2013 : 


leated by the 


Leipski et al. 


20141) . In our templates, these hot structures look more prominent in the Syl than in the 
Sy2. However, in both SEDs, we can still see a peak at 100 /rm, which implies that star 
formation may still be active in these galaxies. The fact that star formation may be present 
in these AGNs is not unexpected, knowing that Sy2s and Syls reside mostly in early-type 
spiral galaxies. 


The question of determining what is the typical level o 
is a difficult one, which is still open to interpretation (e.g.. 


star formation in BLAGNs 


Roseboom et al. 


20131) . In 


QSOs, the problem is even more complicated, because we know very little about the 
morphology of their host galaxies, and, consequently, we do not know what is their SFHs. 
Accordingly, we have chosen two different SEDs from the literature that represent the 
extremes. The SED represents a QSO where the AGN component is predominant over 
the star formation component, while the SED ^6 represent a QSO, also with a strong AGN 
component, but where the level of star formation is comparable to what is observed in 
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ULIRGs, like M82 or Arp 220. Ikeipski et al.l ( 2014 ) determined this SED from a sample of 
QSOs at high redshifts, co 
at high redshifts (see also 


Bertoldi et ah 


2003; 

Wane: et ah 

2008, 

2010 ) 


To determine which SED in our sample of templates better reproduces the W2 — W3 
colors of NELGs with different activity types, we have calculated the colors expected from 
these SEDs at different redshifts. It is well known that due to the increase in redshift, the 
wavelengths that enter a set of filters correspond to a bluer part of the SED. This effect is 
at the basis of the K-correction. In Fig. [TTl the shaded areas show the regions of the SEDs 
that produce the W2, W3, and W4 magnitudes, from z = 0 (limit to the red) to ^ = 0.25 
(limit to the blue). We have verified that due to the different filter responses, the difference 
in width of the filters (especially W3) does not produce any significant difference in colors. 
Therefore, to calculate the different magnitudes, we have used the fluxes at the center 
wavelengths of the hlters at different redshifts (the shaded areas), assuming hlter responses 
of one. To convert the fluxes to magnitudes, we have applied the flux correction factors fo r 
a power law spectrum with spectral index a = 0, which, according to 
is the most general case. 


Wright et ah 


(l2Q10|), 


We show the variations of W2 — W3 colors between z = 0 and z = 0.25 for the SFGs 
in Fig. UHl for the TOs in Fig. [191 aiid for the AGNs in Fig. [201 In all these figures the 
light grey areas correspond to the 5% and 95% color percentiles, whereas the darker areas 
correspond to the 25% and 75% percentiles. Note that the dispersions of the W2 — W3 
colors in all these figures are remarkably small, the darker areas having typical widths 
~ 0.5 —0.7 mag. Therefore, we identify the “typical” SEDs of NELGs with different activity 
types with the SEDs that best reproduce their median W2 — W3 colors at redshift z = 0. 
The K-correction, consequently, would correspond to the color correction implied by this 
SED at different redshifts. 
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For the SFGls and SFG2s, in Fig. [18], the best SEDs at ^ = 0 are the SED #2 
(Arp 220), for a typical starbnrst galaxy, the SED #3 for a typical Sy2 and the SED #6, 
which is the SED of a QSO hosted by a galaxy experiencing a high level of star formation. 
All these SEDs have in common that they imply a high level of star formation in the host 
galaxies, which suggests that this characteristic is a necessary condition to explain the 
position occupied by the SFGs in the MIRDD. 

For the TO Is, in Fig. [191 fhe same three SEDs as for the SFGs fit their median color 
at z = 0, suggesting that their host galaxies are also actively forming stars. However, the 
medians of these SEDs are slightly too high to fit the median of the T02s, and a SED 
giving a value in between SED ^3 and SED ^5 would be better. This would be consistent 
with a lower level of star formation in the T02s as compared to that of the TOls and SEGs. 

Similarly, for the Sy2s and LINERs in Fig. [20] we see a clear shift in median colors 
toward the blue. The best SEDs at z = 0 for the Sy2s are either the SED ^4, which is the 
SED of Mark 231, a Syl galaxy, or SED ^5, which is the SED of a typical QSO with a low 
star formation component. For the LINERs both SEDs, SED ^4 or SED ^5, would only 
fit the reddest of these galaxies, their median being slightly bluer. This is consistent with 
lower star formation in LINERs than Sy2s. 

From the identihcation of the best SEDs at z = 0, we conclude that the different 
W2 — W3 colors of the NELGs with different activity types can be explained by varying 
the contributions of the AGN and star formation components in their SEDs. The AGN 
component is more prominent in the LINERs and Sy2s, whereas the star formation 
component is more important in the TOls and SFGs, while the T02s are between the Sy2s 
and TOls. This result is in good agreement with our analysis of the SFHs of these galaxies. 

But what about the variations of the W2 — W3 colors with the redshift? Examining 
Fig. [TS] Fig. [191 and Fig. [20] the most important characteristics of the NELGs in our sample 
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is that they all show almost constant MIR colors from ^ = 0 up to ^ = 0.25. This explains 
why NELGs with different activity types form specihc regions in the MIRDD. However, 
from the point of view of the SED, we can only reproduce this constancy in colors if it 
includes an AGN component. Obviously, for the Sy2s and LINERs in Fig. [20] this does not 
pose any problem. Gonsequently, adopting SED ^4 the K-correction for the colors of the 
Sy2s at z = 0.25 would be 0.6 mag to the blue, whereas only ~ 0.3 mag would be needed 
adopting SED ^5. Similarly for the LINERs at z = 0.25, a K-correction in color of 0.3 mag 
to the blue would be needed adopting SED #4, while a K-correction of ~ 0.2 mag to the 
red would be needed adopting SED ^5. In both cases, the best £t with lower K-correction 
is SED 7 ^ 5 . However, considering the colors of SED 7 ^ 6 , it seems clear that a QSO SED 
with slightly higher star formation would £t better the Sy2s, while a QSO SED with lower 
star formation would £t the LINERs. 

On the other hand, in Fig. [191 adopting for the TOs the SED #2 for a starburst would 
imply increasingly large K-corrections to the blue for galaxies at z > 0.07, reaching ~ 2.0 
mag at z = 0.25 for theTOls, and slightly less for the T02s. However, this choice of SED 
for the TOs is not unique, since, according to the BPT-VO diagram, TOs are also assumed 
to have an AGN at their center. For example, adopting the SED ^3 would reduce the 
K-correction to the blue to ~ 0.8 mag at z = 0.25 for the TOls, and slightly less for the 
T02s. However, the SED that predicts the lowest K-correction for the TOs is SED ^ 6 , 
producing 0.3 mag to the red for the TOls and 0.5 mag for the T02s. This is the SED 
of a QSO with strong star formation. Therefore, a solution consistent with the average of 
SED 7^5 and SED 7 ^ 6 , would neatly fit the median colors of the TOs at any redshifts, with 
a minimal K-correction up to z = 0.25. But, this would imply increasing the level of star 
formation in TOs at higher redshifts. 


Obviously, the problem of keeping a constant W2 — W3 color as the redshift increases 


28 


is much more critical for the SFGs in Fig. [181 Coincidentally, adopting the SED ^6 for the 
SFGs would predict almost no K-correction up to 2 ; = 0.25 for the SFG2s, and a slightly 
blue correction for the SFGls. But, this solution implies an AGN component is needed in 
the SEDs of these galaxies. Note that this might explain why the colors of the SFGs in the 
MIRDD do not fit those of a black body. Alternatively, a solution based on an average of 
SED 7^1 with SED 7 ^ 2 , with a bias toward low star formation (SED 7 ^ 2 ) at low redshift 
and a bias toward higher star formation (SEDt^I) at higher redshift, may reproduce the 
constant MIR colors of the SFGs without the need for a high K-correction. Adopting this 


solution, star forma 


ion in SFGs would t 


len be expected to increase to the level of ULIRGs 


at higher redshifts (ITorres-Papaqui et al 


2012 


According to our analysis of the SEDs, we conclude that it is possible to reproduce 
the W2 — W3 colors, and the lack of variation with the redshift from z = 0 to z = 0.25 
with minimal K-correct ions, by assuming that: 1) an AGN component dominates the SEDs 
of LINERs and Sy2s, the Sy2s showing a higher level of star formation than the LINERs, 
2) the SEDs of TOs must include both, an AGN and star forming components, with a 
level of star formation decreasing along the sequence T01^T02—)-Sy2, and 3) either the 
SFGs show an AGN with high level of star formation, or star formation in these galaxies 
gradually increases to the level of ULIRGs at higher redshift. In general, these results are in 
good agreement with our analysis of the SFHs, showing that the W2 — W3 colors of NELGs 
with different activity types follow an increase in star formation in their host galaxies. 


^Alternatively, in fiGoziol et al 


20141) it was suggested that a SED with multiple black 


body components could produce MIR colors consistent with a power law. F or the SFGs 


such a SEP could be consistent with a sequen ce of bursts of star formation fiGoziol 


Goziol et al. 


2001 


Torres-Papaqui et al. 


2012 


1996 
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3.4. Comparing BLAGNs with NELGs in the MIRDD 

In Fig. [211 we present the MIRDD for the BLAGNs. We hnd no difference between 
the positions of the Syls and QSOs. They both occupy the same region of the MIRDD, 
which is intermediate between the positions of the LINERs and Sy2s. This is more obvious 
in Fig. [22] where we present the box-whisker plots, comparing the colors of the nearby 
BLAGNs with the median values observed for the NELGs with different activity types. The 
similarities between the nearby QSOs and Syls are conhrmed statistically by the conhdence 
intervals presented in Fig. [23) 1) the Syls and QSOs have similar W2 — W3 colors, whereas 
the QSOs have slightly redder W3 — W4 colors than the Syls, 2) the BLAGNs have redder 
W2 — W3 colors than the LINERs, and are bluer than the Sy2s by about the same amount 
(~ 0.2 mag), and 3) the W3 — W4 colors of the BLAGNs are comparable to those of the 
Sy2s, being much redder than the LINERs. 

We present the MIRDD for the BL Lac objects in Fig. [211 Their positions follow a 
power-law with spectral index a varying between 0 and —1.5, exactly like the LINERs. The 
similarity in MIR colors between the BL Lac objects and LINERs can also be observed in 
Fig. [2H where we show the box-whisker plots for the MIR colors: the BL Lac are bluer 
by ~ 0.4 mag in W2 — W3 than the LINERs (only ~ 0.2 if we consider the BEff using 
w7Gmag), but showing similar W3 — W4 colors. The similarities between the BL Lac and 
LINERs are conhrmed by the conhdence intervals from the max-t tests as presented in 
Fig. [2ni The fact that the LINERs and BL Lac objects have comparable MIR colors is 
important, since these two galaxies constitute the two extremes in terms of AGN luminosity. 

Taken as a whole, one can see that at low redshift {z < 0.25) all the diherent AGNs 
known, namely, the LINERs, the Sy2s and Syls, the QSOs and BL Lac, occupy the 
same region of the MIRDD, their colors tracing a continuous sequence that is consistent 
with the color sequence produced by a power law, Fi, = where the spectral indice, a, 
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varies between 0 and —2, in good agreement with what is seen in the optical-nltraviolet 
part of the spectra of AGNs with different Inminosities. In the MIRDD, one can also 
observe that the W2 — W3 colors of the AGNs become redder along the seqnence 
BL Lac—j-LINER—)-QSO/Syl—;-Sy2. This is different from the AGN Inminosity, which 
decreases along the sequence BL Lac/QSO—)-Syl—)'Sy2—)-LINER. This suggests that the 
sequence in W2 — W3 color is not driven by the AGN component. 


3.5. SED study for the BLAGNs 

In Fig. [271 we show the variations of the W2 — W3 colors of the BLAGNs from z = 0 to 
2 ; = 0.25. For the Syls, there are no clear best hts at 2 ; = 0. The reddest Syls can be htted 
by the SED #2, the SED #3, or the SED #6, exactly like for the SFGs, while the bluest 
one can only be htted by either SED #4 or SED #5, like for the Sy2s and LINERs. On the 
other hand, only SED 7^3 (for a Sy2) or SED would reproduce the colors of Syls at any 
redshift, with minimal K-corrections. These hts seems to imply that the Syl host galaxies 
show a variety of star formation, some being as active as Sy2s, and others much less active, 
similar to the LINERs. This is consistent with their intermediate positions in the MIRDD. 

For the QSOs, the best hts at ^ = 0 are either SED #4 or SED #5. However, SED #5 
would predict smaller K-corrections than SED over the whole range in redshift. The 
diherences in hts suggest that at low redshifts star formation is higher in Syls than in 
QSOs, while it is comparable at high redshifts. The trend for the BLAGNs is thus very 
clear, they both require a strong AGN component, and some level of star formation in their 
host galaxies, with possibly a slightly higher star formation level, on average, in the Syls 
than in the QSOs. 

The case for the BL Lac objects in Fig. [28] is also very clear. None of the six SEDs in 
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our templates sample can reproduce their W2 — W3at2; = 0, these galaxies being much 
too blue in MIR. What we need, it seems, would be a SED with a much higher flux at W2 
than at W3, which, according to our interpretation, would be consistent with a strong AGN 
component in a galaxy with very littl e, or no star formation. To illustrate this point, we 


show in Fig. [29] the SEDs observed by 


Roseboom et ah 


(120131 ). Although the interpretation 


of these SEDs made by these authors do not agree with ours (their models assume that all 
the IR emission is due to the AGN), we note that the SEDs identified as b) and d) would 
have the correct form for the BL Lac objects, showing higher flux at W2 than at W3. 
Indeed, in Fig. [30] one can see that the SED that best reproduces the colors of the BL Lac 
objects is SED d, confirming that for the BL Lac, a strong AGN component is necessary, 
and, obviously (there are no peaks at lOO/rm in these SEDs), no star formation is expected 
in their host galaxies. Again, this is consistent with wh at we know a priori of these objects: 


BL Lac object are 


2000 


Falomoetal 


uminous AGNs in elliptical galaxies flKotilainen et ah 


201 


1998 


Scarpa et ah 


4. Discussion 


According to the BPT-VO diagram, there are only two main sources of ionization of 
the gas in NELGs at low redshift {z < 0.25): young massive stars in star forming regions 
and a SMBH in AGNs. In many AGNs these two sources are known to be operating at 


the same time, which mav complicate their identification (Ho et ah 

.997: Wu et ah 

1998; 

Maoz 1999: 

Vila-Vilaro 

2000 ; 

Kewlev et a 

2001 

Garter et ah 

2001 ; 

Ghen & Zhane 

2006: 

Gonzalez D( 

derado et ah 

2008 

Yuan et ah 

2010 : 

Garoineti et al 

201 ^ 

h Feltre et al. 

2013; 

Goziol et ah 

2014 L So. the question is can we confirm this simple “standard” view in the 


MIR? Our study based on the SFHs and SEDs suggests that the answer is indubitably yes. 
In the MIRDD, NELGs with different activity types, as determined using the BPT-VO 
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diagram, show clear distinctive colors, due to the variation of these two components: a 
SMBH is prominent in the LINERs and Sy2s, whereas star formation increases, as the 
W2 — W3 color becomes redder, along the sequence: LINER—)-Sy2—)-TO—)-SFG. 

Now, when we compare in the MIRDD the BLAGNs with the NELGs, which is 
impossible to do in the BPT-VO diagram, we hnd that the BL Lac objects, the Syls and 
QSOs, fall in the same region of the diagram as the LINERs and Sy2s, their colors forming 
a sequence that is consistent with the colors predicted by a power law, Eh = where 
a varies between 0 and —2, which is standard for AGNs in the optical-ultraviolet. The 
MIRDD, therefore, seems like a robust tool to distinguish between AGNs and star-forming 
galaxies in the nearby universe. In particular, it eliminates the confusion and ambiguities 
about LINERs, clearly showing that they are genuine AGNs. 

Another important characteristic of the MIRDD is the sequence in W2 — W3 colors. 
For the NELGs, including the Sy2s, we have clearly shown that this sequence is due to a 
systematic variation of SFHs with the activity type of the galaxies. Gomparing the BLAGNs 
and BL Lac with the NELGs, we have then found that the color sequence can be extended 
to: BL Lac—)-LINER—)-QSO/Syl—)-Sy2—)-TO—)-SFG. Our analysis of the SEDs suggests, 
therefore, that although QSOs, Syls and BL Lac have a dominant AGN component, 
more comparable to LINERs than Sy2s, their differences in W2 — W3 are consistent with 
different levels of star formation in their host galaxies, which would thus decrease along the 
sequence: BL Lac—)-QSO—)-Syl. This result suggests that the W2 — W3 color trace the level 
of formation in all the emission-line galaxies, irrespective of the importance of the AGN 
component. 

What is the alternative interpretation for the W2 — W3 color sequence? There is a clear 
trend in the literature to favor an interpretation which assumes that in BLAGNs and BL 
Lac objects the MIR emission is dominated by the AGN. This assumption is obviously based 
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on what we see in the optical-nltraviolet part of the spectra. However, according to this 
assnmption, one would naturally expect the MIR colors of these galaxies to follow the AGN 
luminosity, and the sequence should have been then BL Lac/QSO—)-Syl—)-Sy2—)-LINER. 
The problem here is the LINERs. We have shown that, contrary to what was presented in 


the literature (c.f.. 


Wright et al. 


20101) . these AGNs do not experience intense star formation 


activity, as the ULIRGs, and consequently their typical MIR colors, and SEDs are similar 
to those of BL Lac and QSOs, not SFGs, as suggested by the above sequence. The fact that 
we have two types of LINERs, showing similar optical spectra, but different MIR colors due 
to their different levels of star formation, is a clear indication that one cannot apply the 
assumption of a dominant AGN feature in the MIR solely based on their optical-ultraviolet 
spectra. 


The same remark applies to the BL Lac objects. 


are considered to be “pure” AGN (e.g.. 


Plotkin et al 


3asec 


2012 


on their optica. 


Massaro et al. 


spec tra, BL Lac 


we a 

so know that these AGNs are loca 

199§; 

Scarpa et al. 

2000 ; 

Falomo et al. 


ed in non active elliptical galaxies flKotilainen et al. 


20131). However, 


2014J ). which have typical blue colors in the MIR. 


This is exactly what the MIRDD and SED analysis suggests: they have slightly bluer 
W2 — W3 colors than the LINERs in the MIRDD, and their SEDs show that it cannot 
include any star formation components, which is consistent with AGNs in elliptical galaxies. 

Similarly, applying no a priori assumption for the BLAGNs in the MIR, the MIRDD and 
SED analysis suggest that QSOs and Syls are similar kinds of AGNs, where different levels 
of star formation in their host galaxies explain the variations in the MIR. Gonsequently, we 
believe that the most probable interpretation for the W2 — W3 color sequence traced by 
emission-line galaxies in the MIRDD is that star formation increases along the sequence 
BL Lac—»-LINER—»-Q SO/Svl—»-Sv2— »-TO^SFG . This interpretation is consistent with 


previous studies (e.g. 


Hamann fc Ferland 


1992 


McLeod &: Riek 


1994 : 


Brotherton et al. 
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19991: IPaee et al 


2001: Heckman et al.l 2004: Schweitzer et al. 2006: Wild et al. 2007 
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2013 
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brmation of SMBH at the 


Mouri fc Tanigiinhi 
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We conclude that he possibility to compare BLAGNs with NELGs is an important 
advantage of the MIRDD. Such a comparison supports the view that there is a tight 
connection between the a ctivity type of ga l axies in the nearby universe and how these 


es form their stars 

(Dokuchaev 

1991 


kamann & Ferland 

1992 

Bovle & Terlevich 

Franceschini et al. 

I999I; 

Haiman et 

al. 

2004: 

Menanteau 

et al. 

2005; 

Lou & Jiane 


1998 : 


2008 


Haiman 


Treister et al 


et a 


2011 


2007 


Li et ah 


Donoso et ah 


Torres-Papaqui et al. 


2013 


200 ^ 


20121 : 


Bennert et al 


2008 


Torres-Papaqui et al 


Fanid, 

ikis et al. 

2011 


2012 a 1 

n Flovd et al. 

2013 


5. Summary and conclusions 


In 


Goziol et ah 


(120141) ■ a new diagnostic diagram in the MIR was proposed, the 


MIRDD, which is based on WISE data. In the present study we have performed a stellar 
population synthesis analysis on the NELGs and compared their MIR colors with those 
produced by different SED templates. We have also compared the MIRDD of BLAGNs and 
BL Lac objects with the MIRDD for NELGs, and compared their SEDs. The main results 
of our study are the following: 


The MIRDD allows to separate emission-line galaxies according to their activity 


types, by identifying the dominant component in the MIR SED of their host galaxies: 
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an AGN is predominant in the LINERs, the BL Lac objects, the Syls, the Sy2s, and 
in the QSOs, whereas star formation is more important in the TOs and SFGs. 

• The LINERs and BL Lac occnpy the same region of the MIRDD as the LINERs and 
Sy2s, consistent with a power law with spectral index varying between 0 and —2.0, in 
good agreement with photoionization model in the optical-nltraviolet. 

• The W2 — W3 color systematically changes in AGNs, in a way which is independent 
of their Inminosities. The most probable interpretation for this change in color is that 
star formation in the host galaxies of nearby emission-line galaxies increases along the 
seqnence BL Lac^LINER^QSO/Syl^Sy2^TO^SFG. 


Onr resnlts nnderline the importance of star formation in nearby AGNs 
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Fig. 1.— BPT-VO diagram for the NELGs selected for our MIR st udy. The boundaries 


betwe en the NELGs with different activity types are those defined in 


(120131) 


Torres-Papaqui et al. 
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Fig. 2.— Histogram showing the normalized redshift distributions of the different emission 
line galaxies in our sample. 






49 


cn 

D 


E 




cn 

o 


E 


to 



0 12 3 4 5 6 

W2 - W3 [mag] 


Fig. 3.— MIRDD for the Sy2s and LINERs. The distributions correspond to normalized 
density contours around the most probable colors. The line with a positive slope separates 
AGNs from SFGs and the one with a negative slope separates NELGs according to the level 
of star formation in their host galaxies. The upper blue curve is the color sequence traced by 


a power law with different spectra indexes, and t 
a black body at various temperatures (Table 1 in 


le lower red curve is the color sequence for 


Wright et ah 


2010). The color 


K giant is also indicated as an asterisk at the bottom left (also in 


Wright et ah 


or a typical 


20101 ) . 
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Fig. 4.— MIRDD for the TOs. The boundaries and curves are as described in Fig. |3l 
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Fig. 5.— MIRDD for the SFGs. The boundaries and curves are as described in Fig. [3l 
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Fig. 6.— Box-whisker plots of the MIR colors of the NELGs with different activity types. 
The lower side of the box is the 25th percentile, Qi, and the upper side is the 75th percentile, 
Qs- The whiskers correspond to Qi — 1.5 x IQR and Qs + 1-5 x IQR, where IQR is the 
interquartile range. The median is shown as a bar and the mean as a circle. Although barely 
visible, the box-whisker plots include notches, which are V-shape regions, drawn around the 
medians, that have a width proportional to the interquartile range: ±1.58 xIQR/\/N, where 
N is the size of the sample. Samples with overlapping notches have comparable medians. 
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Fig. 7.— Results of the max-t tests in terms of the 95% confidence intervals. This para¬ 
metric test compares simultaneously the differences of mean MIR colors between each pair 
of subsamples (family-wise). A confidence interval including zero indicates no statistically 
significant difference in colors. Negative color differences suggest that the first subsample of 
galaxies in the pair compared are on average bluer than the second subsample of galaxies. 
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Fig. 8.— Variations with redshift of the W2 — W3 colors of galaxies with different activity- 
types and levels of resolntion, according to the WISE photometry flag (ext-flags): 0 (green), 
1 (red) , [2,3,4], (bine), and 5 (yellow). The gaps at z ^ 0.115 and 2 ; ~ 0.148 for the NELGs 


are observational artefacts in SDSS. 
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Fig. 9.— Variations with Petrosian radins of the W2 — W3 colors of NELGs with different 
activity-types and levels of resolntion. The colors have the same meaning as in Fig. |H1 
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Fig. 10.— Comparison of the W3 — W4 colors with the W3G — W4G colors, with different 
WISE photometry flag wAgflg values: 0 (red), 1 (green) , 17 (blue). 
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Fig. 11.— Variations of the SFH in NELGs with different activity types at different time. 
The abscissa gives the age of the stellar populations in the library used for our analysis. The 
continuous curves are the median SFRs and the two dashed curves are the first and third 
quartiles. 








Fig. 12.— Box-whisker plots showing the variations of the two characteristic SFRs in NELGs 
with different activity types: a) SFRyoung, the most recent maximum in star formation rate, 
and b) SFRoid, the maximum star formation rate in the distant past. 
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Fig. 13.— Confidence intervals comparing the two characteristic SFRs in NELGs with 
different activity types 
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Fig. 14.— Relation between the maximum star formation in the distant past, and the age 
of the stars forming this maximum. 
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Fig. 15.— Variations of SFRyoung as a function of the MIR colors of the NELGs with different 
activity types: a) W2 — W3, b) W3 — W4. 
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Fig. 16.— Variations of SFRom as a function of the MIR colors of the NELGs with different 
activity types: a) W2 — W3, b) W3 — W4. 
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Fig. 17.— SED templates of NELGs with different activity types. The SEDs are identihed 
as 7^1 = M82, 7^2 = Arp 220, ^^3 = Sy2, #4 = Syl (Mrk 231), 7)^5 = a QSO with relatively 
mild star formation in its host galaxy, and 7)^6 = Leipski, a QSO with relatively high star 
formation in its host galaxy (see explanations in the text). The shaded areas show the 
regions of the SED that produces the W2 — W3 and W3 — W4 colors in the MIRDD from 
z = 0 (limit to the red) to 2 : = 0.25 (limit to the blue). 
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Fig. 18.— Observed variations of the W2 — W3 colors of the SFGs from z = 0 to z = 0.25. 
The light gray area corresponds to the 5% and 95% percentiles, whereas the darker area 
corresponds to the 25% and 75% percentiles (or hrst and third quartiles). The different 
curves, identified by the numbers in the right margin, correspond to the colors at different 
redshifts, as predicted by the SEDs with the same number in Fig. [T71 
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Fig. 19.— Observed variations of W2 — W3 color of the TOs from z = 0 to z = 0.25. The 
shaded areas and curves are as explained in Fig. [THl 
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Fig. 20.— Observed variations of W2 — W3 color of the Sy2s and LINERs from ^ = 0 to 
z = 0.25. The shaded areas and curves are as explained in Fig. [181 
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Fig. 21.— MIRDD for the nearby BLAGNs. In a) the Syls and in b) the QSOs. The 
boundaries and curves are as described in Fig. 2. 
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95% family-wise confidence level: W3-W4 
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Fig. 23.— Confidence intervals for the comparison of MIR colors of BLAGNs with the colors 
of the Sy2s, the LINERS, the TOls and SFGls: a) W2 - W3 colors, b) W3 - W4 colors. 
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Fig. 24.— MIRDD for the nearby BL Lac objects. The boundaries and curves are as 
described in Fig. 2. 
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Fig. 25.— Box-whisker plots comparing the distributions of MIR colors of the nearby BL Lac 


objects with the median colors of the Sy2s, the LINERS, the TOls and SFGls. 
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95% family-wise confidence level: W2-W3 
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Fig. 26.— Confidence intervals for the comparison of the MIR colors of the nearby BL Lac 
objects with the colors of the Sy2s, the LINERS, the TOls and SFGls. 
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Fig. 27.— Observed variations of W2 — W3 color of the Syls and QSOs from 2 ; = 0 to 
z = 0.25. The shaded areas and curves are as explained in Fig. [THl 
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Fig. 28.— Observed variations of W2 — W3 color of the BL Lac objects from = 0 to 
z = 0.25. The colors correspond to the different level of resolutions according to the WISE 
photometry flag (ext-flags): 0 (green), 1 (red) , 234, (blue), and 5 (yellow). 
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3s fo r QSOs with prominent AGN components, as establsihed by 


(120131 ). The shaded areas show the regions of the SED that produce the 


W2 — W3 and W3 — W4 colors from z = 0 to z = 0.25. 
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Fig. 30.— Comparison of the variation of W2 — W3 color with redshift for the BL Lac using 
the SEDs in Fig. [221 
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Table 1. Results for the SFH of NELGs 


Activity 

Type 

Log(SFRYoung) 
(M0 yr-i) 

Log(SFRoid) 
(Mq yr-i) 

Log(tsFMAx) 

(yrs) 

LINER 

-0.16 

0.66 

8.71 

Sy2 

0.19 

0.64 

9.07 

T02 

0.33 

0.55 

9.20 

TOl 

0.38 

0.44 

9.34 

SFG2 

0.50 

0.28 

9.37 

SFGl 

0.43 

0.23 

9.59 






